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Abstract. A new sufficient condition is proved for the existence of sto- 
chastic semigroups generated by the sum of two unbounded operators. It 
is applied to one-dimensional piecewise deterministic Markov processes, 
where we also discuss the existence of a unique stationary density and 
give sufficient conditions for asymptotic stability. 



1. Introduction 

The development of cell cycle models to account for the statistical prop- 
erties of division dynamics in populations of cells inevitably led to the con- 
sideration of stochastically perturbed dynamical systems [51 fTTl [131 1231 124] . 
These applied considerations have been followed by work on the behaviour of 
Poisson driven dynamical systems in a pure mathematical context p~4l [2Tj . 
More recently other areas of application related to the role of intrinsic (as 
opposed to extrinsic) noise in gene regulatory dynamics P, [T5l [4] have made 
the understanding of stochastic perturbations of dynamical systems of more 
than passing interest. 

We were originally motivated by the work of Lasota et al. [32] who con- 
sidered a (biological) system which produces 'events' and has an internal 
or physiological time in addition to the laboratory time t. We denote this 
internal time by r to distinguish from the time t. When an event appears 
the physiological time r = r e is reset to r = 0. We assume that the rate 
dr/dt depends on the amount of an 'activator' which we denote by a. Thus 
we have 

dr 

(1.1) 37 = V?K V>0. 
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The activator is produced by a dynamics described by the differential equa- 
tion 

da , . 

(1.2) Tt = g (a), 

where g > is a continuous function on an open interval that may or may 
not be bounded. When an event is produced at a time r e and activator level 
a e , then a portion p(a e ) of a e is consumed so the level of the activator after 
the event is then 

(1.3) cr(a e ) := a e - p(a e ). 

We also allow the possibility that the portion p(a e ) depends on an environ- 
mental or external factor so that p is a function of two variables p(a e , 9 e ) 
where 9 e 6 © is distributed according to some probability measure v on 0. 
The solution of (11.21) with the initial condition a(0) = x will be denoted by 

a{t) = ir t x 

and we assume that it is defined for all £ > 0. Then the solution of equation 
(11.1 h with the initial condition r(0) = is given by 

r(t) = / (p(n r x)dr. 



It is reasonable to require that also r(t) is finite for alH > 0. 

Lasota et al. [13] studied the statistical behavior of a sequence of such 
events occurring at random times 

= t < h < ■ ■ ■ < t n < ■ ■ ■ 

and denoted a n = a(t n ) to find 

ftn 

(1.4) a n+ i = T(a n ,r n ), where r n = ^>{'Ks-t n - 1 an-i)ds 

were exponentially distributed independent random variables, giving a rela- 
tion between successive activator levels at event occurrence and studying a 
discrete time system with stochastic perturbations by the r n . Here we have 
extended these considerations to a continuous time situation by examining 
what happens at all times t and not merely what happens at to, t\, t2, ■ ■ ■ ■ 
Thus we arrive at a continuous time piecewise deterministic Markov pro- 
cess, whose sample paths between the jump times to,ti,--- are given by the 
solution of (II .21) and at the jump times the state of the process is selected 
according to a jump stochastic kernel, which is the transition probability 
function for (jl .31) . This leads us to study evolution equations of the form 

Ou d 

(1.5) — = A u — ipu + P(<pu), where A u(x) = — — (g(x)u(x)), 

CJL (JjJb 
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on the space of integrable functions L 1 , where P is a stochastic operator on 
L 1 corresponding to the jump stochastic kernel and cp need not be bounded. 
We supplement f 1 1 . 5 H with the initial condition u(Q) = uq which is the density 
of the distribution of the initial amount of the activator. 
Let us write 

(1.6) Au = A u — ipu and Cu = Au + P(<pu). 

If <p is unbounded then C is the sum of two unbounded operators, so the 
existence and uniqueness of solutions to the Cauchy problem in L 1 is prob- 
lematic; (11.51) may have multiple solutions [2]. We make use of perturba- 
tion results for positive semigroups on L 1 -spaces which go back to [10] (see 
Section [3]) , from which it follows that the operator C has an extension C 
generating a positive contraction semigroup {P(t)} t >o provided that the op- 
erator A is the infinitesimal generator of a positive contraction semigroup 
on L 1 and C is defined on the domain of A. In general, if the closure of 
C is the generator C of {P{t)} t >o then the Cauchy problem is uniquely 
solved and {P(t)} t >o is a stochastic semigroup. In Section [3] we prove a 
new sufficient condition for uniqueness and in Section [5] we show that if the 
discrete process has a strictly positive stationary density then uniqueness 
holds. This simplifies the analysis of (11.51) when compared with the ap- 
proach in [2], and allows us to investigate both the uniqueness of solutions 
and their asymptotic properties. 

The outline of this paper is as follows. We recall basic definitions and fun- 
damental theorems from the theory of stochastic operators and semigroups 
in Section [2] and perturbation results for positive semigroups on L 1 -spaces 
in Section El which closes with the proof of our main general result (The- 
orem [3]) . In Section [4] we prove that the operators A and A defined on 
suitable domains are generators. In Section [5] we show the applicability of 
Theorem [3] to equation (11.51) when P is an arbitrary stochastic operator and 
also give sufficient conditions for asymptotic stability. In Section [6] we let 
the operator P have a definite form that fits directly into our framework, 
and give several concrete examples drawn from work on the regulation of 
the cell cycle as well as classical integro-differential equations. In Section [7] 
we extend our results to the situation in which there is degradation (as 
opposed to growth) and illustrate their applicability using models for the 
stochastic regulation of gene expression. 

In a companion paper [22], these and other results are placed in the gen- 
eral context of semigroup theory and probability theory with applications to 
piecewise deterministic Markov process without 'active boundaries'. There 
we also prove that when the semigroup {P(t)} t >o is stochastic then (1 1.5ft is 
the corresponding evolution equation for densities of such processes. 
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2. Stochastic operators and semigroups 

Let (E,£,m) be a a-finite measure space. We denote by D the set of all 
densities on E, i.e. 

D = {u £ L 1 : u > 0, |H| = 1}, 

where || • || is the norm in L 1 = ^(E, £,m). A linear operator P: L 1 — > L 1 
such that -P(-D) C D is called stochastic or Markov [12] . 

Let P: L 1 — > L 1 be a stochastic operator. A density u is said to be 
invariant or stationary for P if Pw = w. We say that P overlaps supports 
if for every u,v £ D there is a positive integer n > 1 such that 

m(suppP n u n suppP n v) > 0, 

where the support of u £ L 1 is defined up to a set of measure zero by the 
formula suppw = {x £ E : u(x) ^ 0}. Note that if P overlaps supports then 
it can have at most one invariant density [HI see the proof of Corollary 1]. 

Let J: E x £ — > [0, 1] be a stochastic transition kernel, i.e. J~(x,-) is 
a probability measure for each x £ E and the function x \— > J{x, B) is 
measurable for each B £ £, and let P be a stochastic operator on L 1 . If 

(2.1) / J(x,B)u(x)m(dx) — / Pu(x)m(dx) for all B £ £, it £ D(m), 
Je Jb 

then P is called the transition operator corresponding to J . 

If J(x,B) = l T -i( fl )(x) for x £ P, P £ f, where T: P -> P is a 
nonsingular measurable transformation, i.e. m(T _1 (P)) = for all P £ £ 
such that m(B) = 0, then there exists a unique stochastic operator P on 
L 1 satisfying (|2.ip and P defined by (12. ip is called the Frobenius- Perron 
operator corresponding to T. 

A stochastic operator P on L 1 is called partially integral or partially kernel 
if there exists a measurable function p: P x P — > [0, oo) such that 

p(x,y) m(dy) m(dx) > and Pu(x) > / p(x,y)u(y)m(dy) 
' e J e Je 

for every density w. If, additionally, 

/ p(x, y) m(dx) — 1, y £ P, 

then P corresponds to the stochastic kernel 

J(x,B)= / p(y,x)m(dy), xeE,Be£ 
Jb 

and we simply say that P has kernel p. 

A strongly continuous semigroup {P(t)} t >o on L 1 is called a stochastic 
semigroup or Markov semigroup if P(£) is a stochastic operator for all t > 0. 
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A density is called invariant or stationary for {P(t)} t > Q if P(£)w* = 
for every t > 0. 

A stochastic semigroup {P(t)} t >o is called asymptotically stable if there 
is a stationary density u* such that 

lim \\P(t)u-u4 = for u e D 

t— >oo 

and it is called partially integral if, for some to > 0, the operator -P(io) is 
partially integral. 

Theorem 1 ([12]). Let {P(t)} t >o be a partially integral stochastic semi- 
group. Assume that the semigroup {P(t)} t > has only one invariant density 
it*. If u* > a.e. then the semigroup {P(t)} t >o is asymptotically stable. 

3. Perturbation results in L 1 

Let (A ,V(A )) be the infinitesimal generator of a stochastic semigroup, 
<f > be a measurable function, and 

L^p = {u E L 1 : / ip(x)\u(x)\m(dx) < oo}. 
Je 

Let P be a stochastic operator on L 1 and let the operators A and C, as 
given in (11.61) . be defined on T>(A) C P(^4 ) HL^. Assume that the operator 
(A, P(A)) is the infinitesimal generator of a positive strongly continuous 
contraction semigroup {S(t)} t > on L 1 . Then it is known [TQl EH], [TJ, [2] that 
there exists a positive strongly continuous contraction semigroup {P(t)} t >o 
on L 1 satisfying the following: 

(1) the infinitesimal generator C of {P(t)} t > is an extension of the 
operator C, i.e. V(A) C V(C) and Cn = Cn for u G X>(A); 

(2) if {P(t)} t >o is another semigroup generated by an extension of C 
then P{t)u > P(t)u for all u e L 1 , u > 0, i.e. {P{t)} t > is the 
minimal semigroup; 

(3) the generator C is characterized by 

N 

(3.1) R(X,C)u— lim R(X,A)y^(P((fR(X,A))) n u, ueL\\>0, 

71=0 

where R(X, •) is the resolvent operator; 

(4) the semigroup {P(t)} t >o satisfies the equation 

(3.2) P{t)u = S{t)u+ [ P{t- s)P{ V S(s)u)ds, ueV(A). 

Jo 

We can not conclude, in general, that the semigroup {P(t)} t > is sto- 
chastic pIU Example 4.3]. Discussing various conditions for this to hold has 
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been a major objective of study pm [25], [U d [2] and leads to the following 
result. 

Theorem 2. If for some A > 

(3.3) lim || (P(ipR(\, A))) n u\\ = for all u G L 1 

n— »oo 

then {P(t)} t >o is a stochastic semigroup and its generator C is the closure 
of the operator (C,T>(A)). 

We now prove our main result. 

Theorem 3. If for some A > there is v G L 1 such that v > a.e. and 
P(<pR(\, A))v < v, then condition (13.31) holds. 

Proof. We can assume that v is a density. Set K\ = P(tpR(X,A)). The 
operator K\ is a positive contraction. Since K\v < v, the sequence K™v is 
strongly convergent in L 1 to a fixed point u* of the operator K\. From (2J 
Theorem 4.3] it follows that = 0. For any density u we have 

< K\Uk < kK™v, where Uk = min{u, kv}, k > 1, n > 1. 

Since ||ttfc — it|| — >• as — > oo, this completes the proof of fj3. 3f) . □ 

4. The semigroup {£(£)}i>o 

Let E be an open interval in M, bounded or unbounded, £ = 13(E) and 
m be the Lebesgue measure. We shall denote by L\ oc the space of all Borel 
measurable functions on E which are integrable on compact subsets of E 
and by AC the space of all absolutely continuous functions on E. We assume 
from now on that E = (do,di), where — oo < do < d\ < oo, g: E — > R is 
a continuous strictly positive function, and if G L\ oc is nonnegative. In this 
section we study the first order differential operators A and A which are 
meaningful for any function u G L\ oc for which gu G AC. We will define the 
operators A and A on suitable domains T>(A ) and T>(A) so that they are 
generators of corresponding semigroups as described in Section [3l 

Since l/g,ip/g G L\ oc , we can define 

(4.1) G(x) = [ ^—dz and Q(x) = [ ^-dz, 

Jxo 9\z) J X1 g(z) 

where x = d and X\ = d when the integrals exist for all x and, oth- 
erwise, Xq, X\ are any points in E. The function G is strictly monotonic, 
continuously differentiable on E, G(d Q ) G {0, — oo}, and G(d\) is either fi- 
nite or equal to +oo. The function Q is monotonic with Q(d ) G {0, — oo} 
and Q(di) is either finite or equal to +oo. G is invertible with G~ x well 



DYNAMICS AND DENSITY EVOLUTION IN GROWTH PROCESSES 



7 



defined on G(E). If G(E) ^ R, then we extend G 1 continuously so that 
G-\R \ G(E)) = {d , rfi}. The formula 

a(t, x) = G~ 1 (G(x) + t), xeE,teR, 

defines a monotone continuous function in each variable with values in 
[—00,00] which is a solution of (11.21) . If G(di) = +00 then n t x = a(t,x) 
is finite for all t > 0, x G E, and n t (E) C E, t > 0. In the case when 
|G(d )| = G(d\) = 00 the value a(t,x) is finite for all t e 1 and x G E, so 
that we have, in fact, a flow n t on E such that 7r t (E) = E. 
For t > we define the operators Po(t) and S'(i) on L 1 by 

(4.2) P (tW x ) = l E (7r_ t x)u(7r- t x) 9 ^~f\ x <E E,u <E L 1 
and 

(4.3) S'(t)u(x) = e Q{w - tx) - Q{x) P {t)u{x), x£E,ueL\ 

Theorem 4. IfG(di) = +00 t/ien {-Po(^)}t>o is a stochastic semigroup and 
{S(t)} t >o is a positive strongly continuous contraction semigroup on L 1 . 

Proof. Po(t) is a stochastic operator because it is the Frobenius-Perron oper- 
ator for the transformation x 1— >■ ir t x. Since Q is nondecreasing and t 1— >■ 7r t x 
is increasing, we always have 

e^-^-^le^x) < l E (7r_ t x). 

Hence S(t) is a positive contraction. To check the semigroup property 
observe that if x G E and n- s - t x G E, then 7r s (7r_ s _ t x) G £7, by assumption, 
and thus 7r s (7r„ s „ t x) = 7r„ t x G -E. Furthermore, if x G and 7r_ 4 x G i? then 
7T-s(7T-^) = T^-s-tX- Hence 

l E {-K-tx)l E {^- s {^-tx)) = l E {^- s -tx), 

which shows that S(t)S(s)u(x) = S(s + t)u(x) for t, s > 0, x G E, u G L 1 . 
Finally we must show that {S(t)} t >o is strongly continuous. Let u G C C (E), 
where C C (E) is the space of continuous functions which are equal to zero 
near boundaries. For every x G E and all sufficiently small t > we have 
l^(7r_ 4 x) = 1 and 7r_ 4 x — > x as t j 0. Consequently, 

lim SYtWx) = w(x) x G E, 

which, by the Lebesgue dominated convergence theorem, implies 

lim\\S(t)u-u\\ = 0. 

Since the set C C (E) is a dense subset of L 1 , this shows the strong continuity 
of the semigroup {S(t)} t > . If we take ip = then S"(t) = Po(0 an d we 
recover the claim for {Po(t)} t > . □ 
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Now we identify the generators of the semigroups from Theorem 0J The 
maximal domain of A in L 1 consists of all functions u G L 1 such that 
A u G L 1 . Then the integrability of A u implies existence of the finite 
limits 

(4.4) lo(u) := lim g(x)u(x) and h(u) := lim g(x)u(x). 

xid x^di 

A necessary condition for A to generate a stochastic semigroup is that the 
limits U are equal. Recall that u G L 1 ^ if and only if u G L l and (pu G L 1 . 

Theorem 5. If G(di) = +oo then the operator A Q defined on the domain 

(4.5) V(A ) = {u <E L 1 : gu <E AC, A u G L 1 , lim g(x)u(x) = 0} 

xldo 

is the generator of the semigroup {-Po(^)}t>o an d the operator A defined on 
V(A) = V(A Q ) n L]p is the generator of {S(t)} t > . 

Before we give the proof of Theorem [o] we first provide general properties 
of the operator A. Let A > 0. Define the function r\ : E x E — > [0, oo) by 

e Qx{y)-Q\(x) 

(4.6) r x (x,y) = l {do , x) (y) — , where Q\(z) = XG(z) + Q(z), 

g{x) 

and the positive linear operator R\ \ L 1 — > L 1 by 



(4.7) R x v(x) = r x (x,y)v(y)dy. 

J d 

Lemma 1. Let A > 0. The operator R\ satisfies 

\\\Rxv\\ < \\v\\, veL 1 . 

For every v G L 1 we have gR\V G AC and the function u = R\V is a 
particular solution in L 1 of the equation 

(4.8) Am - Au = v. 

Proof. R\ is an integral operator with nonnegative measurable kernel r x . 
Observe that 

A 

' g(x)' 

Thus for every y G E we have 



Ar A (x, y) < l {<khx) {y) — eX{G(v) ~ G{x)) for x, y G E. 



A / r x {x, y) dx < 1 - e -KG{d x )-G{ y )) < ^ 



because G(di) > G(y), showing A||-Ra|| < 1- Since e Q G AC, we can write 

Au(x) = -e~ Q{x) ^-(g(x)u(x)e Q(x) ), 
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for all u such that gu G AC. We have 

g(x)R x v(x) = e~ Qx{x) / e Qx{y) v(y) dy 

Jdo 

and Q\(y) < Q\{oi) for every y < a < d\, thus the function e Qx v is inte- 
grable on intervals (d ,a\ for every a < d\. Hence, gR\v G AC and R\v 
satisfies {4l§. □ 

Lemma 2. If Q\{d x ) = +00 then RxiL 1 ) C V(A) and 

R x (\u-Au) = u for ueV(A). 

Proof. First we show that R\{L l ) C X?(A). Let t> G L 1 . Then ^i? A w G AC, 

\g(x)R\v(x)\ < / \g{x)r x (x,y)\\v(y)\ dy for x E E, 

Jdo 

and |(7(x)r A (x, ?/)| < l( do for all x,y E E. From the definition of r A and 
the assumption Q\(di) = +00 it follows that 

lim \g(x)r\(x, y) \ = for y E E. 

x—>di 

By the Lebesgue dominated convergence theorem, li(R\v) = 0, i — 0, 1, for 
all f G L 1 . Since |t>| G L 1 , the function -RaM is a particular solution of 

(\ + <p)R x \v\ = \v\ + A (R x \v\) 

and li(R\\v\) = for i = 0, 1. Hence 

(A + (f(x))R\\v\(x) dx = / \v(x)\dx, 

do J do 

which shows that ||</?-RaMII < 00 and R\v G L 1 ^. Finally, from (14.81) it 
follows that 

Aq(R\v) = (A + ^)i? A u — v £ L 1 . 
Now let m G £>(A) and v := Xu — Au. Since and A u G L 1 , we 

have v E L 1 and i?Af G ^(^4). Hence w := u — R x v G £>(A) and Au; = Aw. 
The general solution w of this equation is of the form 

e -Qx(x) 

w(x) = c — —— , x G E, 

where c is a constant. Thus 

g(x)(u(x) - R\v(x)) = ce~ Qx ^ x) for x G E. 

If |<5(<io)| = +00 and c 7^ 0, then ce _QA (A + y)/^ is not integrable, which 
contradicts u—R\v G £>(A) and gives c = 0. If Qa(^o) = then (u—R\v) = 
c, which also gives c = and shows that u = R\v. □ 
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Proof of Theorem^ Since G(di) = +00, the assumptions of Lemma [2] hold 
even when <p = 0. Thus it is enough to show that (A, T>(A)) is the generator 
of {S(t)} t >Q. Observe that L 1 = Im(A — A), by Lemma El Since (14.81) has 
a unique solution u = R\v, we have R\ = (A — A)" 1 for A > and 

|| Aw - Au\\ = \\v\\ > \\\R\v\\ = X\\u\\ for u e V(A). 

The operator \R\ is a positive contraction. By the Hille-Yosida theorem, 
the operator A is a generator of a positive contraction semigroup. Let A be 
the generator of the semigroup {S(t)} t >o- It remains to prove that A = A. 
First, we will show that 

R{\, A)v = R x v for v G C C (E), v > 0. 

Let v G C C {E), v > 0. We have 

R(X,A)v = lim / e~ Xs S{s)vds, 

where the integral is an element of L 1 such that 



-A.s 



S(s)vds ] (x) 



-A.s 



S(s)v(x) ds. 



Let x G E and t > 0. Define s*(x) = sup{s > : n_ s x G E 1 } and note 
that 7r_ s »(a;)X = do- Making use of the formula for S(t)v and the fact that 
G{ji^ s x) = G(x) — s when 7r_ s x G -E leads to 



-As 



S(s)v(x) ds 



tAs,(x) 



&-.x)-Q x {x) v{^- s x)g{n_ s x) 



ds. 



By a change of variables, we obtain 



-A.s 



S(s)v(x) ds — R\v(x) 



< 



1 



w(t, x) 



where 



w(t, x) 



TT-tX 



e* (2) w(z) dz, if t < s m (x); 

0, ift>s*(x). 

We have w(t, x) [ as t ] 00, thus 

-t 



lim 

t— >oo 



e- Xs S(s)vds - R x v 



0. 



Since C C (E) is a dense subset of L 1 and both operators XR(X, A) and XR\ 
are positive contractions, they are identical. This shows that T>(A) = T>(A) 
and (A — A)u = (A — A)u for u G ^(^4), which completes the proof. □ 

Remark 1. Observe that if G(do) = —00 then the domain of A is 
V(A Q ) = {ueL 1 :gue AC, A u G L 1 }. 
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5. Asymptotic properties 

In this section we assume that g > and G(dx) = +oo. Let P: L 1 — > L 1 
be a stochastic operator and let C be the operator 

Cu(x) = — — (g(x)u(x)) — ip(x)u(x) + P((pu)(x), x G (do, di), 

defined on V(A) = V(A ) n h\ with V(A ) as in (fHHh By Section El it 
follows from Theorems 0] and [5] that there is a positive, strongly continuous, 
contraction semigroup {P(t)} t >o on L 1 whose generator is an extension of 
the operator (C,V(A)). In this section we give sufficient conditions for 
{P(t)} t >o to be a stochastic semigroup and study its asymptotic properties. 
Define the operator i? on V(Rq) = {v G L 1 : R v G L 1 } by 

Rov(x) = r (x,y)v(y)dy J where r (x, y) = l {do>x) (y)— e Q{y) ~ Q{x) . 

In general, Ro may be an unbounded operator. We have 

Mi 

ip(x)R v(x) = / (f(x)r (x,y)v(y)dy 

Jdo 

for every t> G V(Ro). Since Q is nondecreasing, 

^ ^(x)ro(z, y) rfx = ^l e Q ^- Q ^dx = 1 - e 0(w)-0(*) < i 

d Jy 9{X) 

for every y E E. Thus the operator ipR can be uniquely extended, with 
the same formula, to a positive contraction on L 1 . Observe that ipRo is 
stochastic if and only if Q(d\) = +oo. 
Define the operator K by 

(5.1) Ku = P((pR )u for u G L 1 . 

Since P is a stochastic operator, K is stochastic if and only if Q(di) = +oo. 

Theorem 6. Assume that Q(di) = +oo. 

// the operator K has an invariant density v * > a. e. then {P(t)} t > 
is a stochastic semigroup. Moreover, if {P(t)} t > is partially integral and 
RqIU G L 1 then {P(t)} t > is asymptotically stable with invariant density 
u* = RoV m /\\RoV*\\. 

Conversely, if {P(t)} t >o has an invariant density G T>(A) then the 
density P(<pu*)/\\(pu*\\ is invariant for the operator K. 

Proof. The operator Ri = R(1,A) as defined in (14 .7K satisfies (<pRi)u < 
(<pRo)u for every density u, which is due to the fact that ip > and r\ < r . 
Since P is a positive operator, we have P(<pR(l, A))v* < v* and {P(t)} t > is 
a stochastic semigroup, by Theorems [3] and El Now let C be the generator 
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of {P(t)}(> . From ( 13. II) it follows that R(l, C)u > Riu for every density it. 
Since there is a b(u) G -E such that Riu(x) > for x > b(u), the stochastic 
operator R(l, C)u overlaps supports. Recall that w* is a stationary density 
for the semigroup {P(t)} t > if and only if w* G T>(C) and Cm* = 0. We 
have 



It is easily seen that u* G £>(A) and An* = — u*/||i2oU*||. Hence Cu* = 
and Theorem [T] applies. 

Finally, suppose that P(t)w* = u* for all £ > with u* G P(-A). Since 
T>(A) C and Cu* = Au* + P(<pu*), we obtain —Au* = P(<pu*). Thus 
w := P(<pu*) G L 1 and Kt> = P((pRov) = —P(<pRoAu*). It is easily seen 
that RoAu* = — u* and we recover the claim. □ 

Remark 2. Observe that when <p{x) > b > for all x, then the operator 
Rq is bounded, thus Rq = R(0, A) = —A~ l and RqV* is integrable. More- 
over, if ip is a constant function, cp = b, then (C,T>(A)) is the generator 
of {P(t)} t >o, by the Phillips perturbation theorem, and K = P(bR(b, A Q )) . 
In that case, the relation between the invariant densities for the operator 
K and the semigroup {P(t)} t >o are a consequence of the equation Cu* = 
which is now A u* — fat* + bPu* = (see also [T4] ). 

General sufficient conditions for existence of invariant densities for sto- 
chastic operators have been summarized by [321 Section 5] and [20]. We now 
discuss when a stochastic semigroup {P(t)} t >o is partially integral. Let P 
be the transition operator corresponding to a stochastic kernel J . If there 
is a Borel measurable function p: E x E — > [0, oo) such that 



then {P(t)} t >o is partially integral pU [H]. Now, if P is the Frobenius- 
Perron operator corresponding to a nonsingular transformation a: E — > E 
then J~(x,B) = 1b(ct(x)), x G E, B G B(E), and we have the following 
result. 

Proposition 1 ([El US]). Let a: E — > E be continuously differentiate with 
a'(x) 7^ for almost every x G E and P be the Frobenius- Perron operator 
corresponding to a. Assume that the semigroup {P(t)} t > is stochastic. If 
there is x G E such that ip is continuous at x, <f{x) > 0, and g(a(x)) ^ 
o~'(x)g(x) then {P(t)}t>o is partially integral. 



v* = P(<p>R )v* 



\\RovJ\P((pu*). 




p(x,y)ip(y)dydx > and J{x,B)> / p(y,x)dy, B G B(E), 




6. Specific examples 



As mentioned in the introduction, we were originally motivated by the 
work of Lasota et al. |13| who, in turn, were trying to understand the ex- 
perimentally observed asymptotic properties of cell property densities in 
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cellular populations in a growth phase. This led, naturally, to a consider- 
ation of dynamics such that g > for x G (do, di) as we have considered 
in Sections [4] and [5J and there are a number of concrete situations in this 
category to which we can apply our results. This section illustrates some of 
these. Several of the examples are drawn from the field of cell cycle kinetics, 
while others illustrate the necessity of certain assumptions. 

Suppose first that the reset function is given by a(x) = x — p(x), where p 
is a continuously differentiable function with p'(x) < 1. Then the transition 
operator P is of the form Pu(x) = X'(x)u(\(x))1e(\(x)), where X(x) = 
is the inverse of a, and we have the evolution equation 

du(^x)_ = J(g(x)u(t,x)) _ y(x)M(t>a . )+M(tiA(a . ))y(A(a . )) y (a . )lg(A(a . )) . 

If Q(d\) = oo then the stochastic operator K as defined in (15. ID is of the 
form 

(6.1) Ku{x) = \\x)^^l (do4l) (\(x)) [ X(X) e^'^Mv) dy, 

9\ A \ X !) Jdo 

which is the transition operator corresponding to (|1.4p . where T(a, r) = 
a{Q- l {Q{a)+r)). 



Example 1. Lasota and Mackey [TT] considered a very general cell cycle 
model for the evolution of the distribution of 'mitogen' at cell birth in which 
g was a C 1 function on [0, 21), such that g(x) > forx > and G(2l) = oo, 
where I is finite or not. Further (f was a continuous function on [0, 21) such 
that <p(0) = 

liminf tp(x) > when I < oo and liminf q(x) > 0, 

x—>2l x— >oo 

where q(x) = ip(x)/g(x). In their model X(x) = 2x. They were able to 
show that successive generations had densities evolving under the action of 
a stochastic operator 



Ku{x) = 2q{2x) 



and that K is asymptotically stable. 



p2x 


p2x 


I exp 


— / q(z)dz 


to 


- Jy 



u{y) dy for < x < I 



Example 2. Building on this model Mackey et al. [IE] took g(x) = x{2 — 
x)/b and <p(x) = S(x — l)l(i ]2 )(x) with b, S > and x G (0,2) to fit a 
number of in vitro cell cycle data sets. With these choices for g and ip it is 
straightforward to show that the unique stationary density of mitogen is 

v*(x) = SV2 5b (a;-i)[x(i-2)] (5b/2) - 1 for xG(|,l). 
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Assume that E = (d , oo) and Q(oo) = oo. We can rewrite the operator 
K given by ( 16.11) in the form 

(6.2) Kuix) = [ X(X) ~ ( e Q ^- Q ^ x ») u(y) dy. 

Jdo ox 

Asymptotic properties of this operator have been well studied [81 [TBI l3l [15]. 
Proposition 2 ([8j). Assume that ip(x) > for x > d and Q(d ) = 0. // 

(6.3) liminf (Q(X(x)) - Q(x)) > 1 

x— >oo 

then K as defined in ( 16.21) has a strictly positive invariant density and if 
Q(\(x)) — Q(x) < 1 for all x > do it has no invariant density. 

Remark 3. The assumption Q(d ) = can not be omitted in Proposition^ 
as the following example shows. 

Let E = (0, oo), A(x) = 2x, and Q(x) = blogx. The operator K is now 

b ^ 



If we take u(x) = 1/x then Ku(x) = u(x) for all b > 0. This shows that 
K has a subinvariant function which is strictly positive and not integrable. 
Since K overlaps supports, it has no invariant density [T8l Remark 6] for 
any b > 0, but condition (16.31) holds whenever b log 2 > 1. 

Example 3. Consider the following functions 

g(x) = k, k > 0, f(x) = px a , p > 0, a > —1, X(x) = 2x, x > 0. 

We have G(x) = kx and Q(x) = bx a+1 /(a + 1), where b = p/k. The 
operator K has a strictly positive invariant density, by Proposition® Thus 
the semigroup {P(t)} t > is stochastic, by Theorem 0, and it is partially 
integral, by Proposition^ since g(x/2) ^ g(x)/2 for all x. The domain of 
the operator A is 

POO 

V(A) = {u e ^[0,00) : u(0) = 0,ue AC,u G L 1 , / x a \u{x)\dx < 00}. 

Jo 

The stationary density w* G T>(A) for the semigroup {P{t)} t >o is a solution 
of the equation 

u\x) = -bx a u(x) + 2b(2x) a u(2x), x > 0, 

and is given by 

00 2 a+1 
= ^2c n e~ Q{2 " x \ where c n = _ 2n(a+1) c ra -i, n > 1, 

n=0 

and Co is a normalizing constant. The stationary density v* for the operator 
K is given by v*{x) = cx a u*(2x), where c is a normalizing constant. 
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We continue with the above example, but now we take a = — 1 and show 
that for certain values of the parameter b the semigroup is stochastic and 
for others it is not. Observe that we have Q(x) = blogx, thus the operator 
K is the operator from Remark [3] and it is not asymptotically stable. 

Example 4. Let the functions g and A be as in Example Let Q{x) = 
blogx, where b = p/k. We have 

u -2kx r2x 

P{yR x )u{x) = -— I / e k yy b u(y)dy. 
z x Jo 

If we take u(x) = x^~ x e~ kx then u G L 1 for [3 > and 

b2P 

PfrRiWz) = VT ^u(x)e- kx . 

Assume thatb\og2 < 1. Since we can find (3>0 such that b2? < b+(3, the 
operator P(ipRi) has a subinvariant strictly positive density, which shows 
that the semigroup {P(t)} t > is stochastic, by Theorem^ 

Assume now that b log 2 > 1 and take k — 1. If we go back to (|1.4p then 



in = \a n -^' b and HO = 



2 v nj 2»t V6-7 

We can find 7 < 1 such that 62~ 7 < b — 7. Since t n = a n — a + Y17=i a i> 
this shows that sup n E(t^) < oo, so that the process is defined only up to a 
finite random time and {P(t)} t >o can not be stochastic. 



Example 5. Tyson and Hannsgen [24| in their cell cycle model consider a 
special case of the model of [11] in which they let E = (a, 00), where a < 1 
and consider the following functions 

g(x) = kx, <p(x) = I °' X x > \ > a ( x ) = (JX - 

They show that the unique steady state is given by 

r — 1 fx: 

v*{x) = - 

a \o~ 

where the exponent r > 1 must satisfy 

1 J) 

b — (r — 1) = bcr'" 1 and bin — > 1, where b = — . 

a k 

In this example we have g(o~x) = o~g(x) for allx and the semigroup {P(t)} t >o 
is not partially integral. Although it has a unique strictly positive stationary 
density, it is not asymptotically stable due to a possible synchronization [5]. 
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We conclude this section with an example when the reset function a 
depends additionally on an external variable. Let E = (0, oo), 6 = (0,1), 
v be a measure on 6 with a density h, and the reset function a be of the 
form o~(x, 6) = x — 6x. Then the transition operator P has the kernel p 

p{x,y) = l( 0tX ){y)4}(-)-, where ^(0) = h(l - 9), 

\y J y 



and the evolution equation is 



du(t,x) dg{x)u{t,x) f 00 /x\tp(y) 

(p{x)u{t,x)+ ip[-\ u(t,y)dy. 



dt dx J x \y 

The operator K has the kernel 

(6.4) k(x,y)= ij(-)^-e Q{y) - Q{z) dz, x, y G (0, oo). 

Jmax{x,y} ^ z ' z 9\ z ) 

Example 6. Suppose that <p(x)/g(x) = bx a for all x > 0, where b > 0, 
a > —1. We have Q(x) = bx a+1 /(a + 1). We provide the form of the 
invariant density for K when if}(z) = (3z^~ l for z G [0, 1] and (3 > 0. It is 
easily seen that the invariant density for the operator K is of the form 

W _/9-i„-Q(*) - P 



v (x) = - x^ _1 e" yw t 

RqV*(x 



(a + l) 1+ Tr( 7 ) ' ' a + l 

We have 

xv*(x) 



(3g(x) • 

If RqV* G L 1 then {P(t)} is an asymptotically stable stochastic semigroup 
by Theorem^ For example when a = then v* is the gamma distribution, 
while if a = (3 = 1 then 

= — bx 1 . 

7. Decay instead of growth 

We assumed in Sections [4] and [5] that g(x) > for x G (d ,di), and 
illustrated the applicability of our results to concrete situations in Section 
[61 In this section we discuss a situation when instead of growth there is 
degradation, so now we suppose that g(x) < for x G (d ,di). Results 
analogous to those of Sections [4] and [5] with similar proofs are valid in this 
case, and we illustrate the applicability of these to models for the stochastic 
regulation of gene expression. 

Let the functions G, Q be defined as in (14.11) . Observe that now G is 
decreasing and Q is nonincreasing. If G(d ) = +oo then n t (E) C E for all 
t > and 1 1 — > ir t x is decreasing. Let {Po{t)}t>o and {S(t)} t >o be as in ( 14.21) 
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and (14.31) . Thus, if g < and G(d ) = +oo the conclusions of Theorem [4] 
remains valid. The analogue of Theorem [5] with the same method of proof 
reads as follows. 

Theorem 7. If g < and G(do) = +oo then the operator A defined on 
the domain 

(7.1) V(A ) = {u <E L 1 : gu e AC, A u G L 1 , lim g{x)u{x) = 0} 

xldi 

is the generator of the semigroup {Po(t)} t >o and the operator A defined on 
V(A) = V{Aq) n L\ is the generator of {S(t)} t > . 

Moreover, the resolvent of the operator A is of the form 

R(X, A)v(x) = [ 1 —ji—e Q ^- Q ^v(y) dy, v G L 1 . 

Jx \g{x)\ 

The operator R is now defined with the help of the function 

r (x,y) = l {Xidl) (y)-±-e^-^ 
\g(x)\ 

and the assertions of Theorem [6] remain valid under the assumption that 
g < and G(d ) = Q(d ) = +oo. 

In particular, if the operator P has kernel p then the operator K has the 
kernel k 

k(x,y)= [ V p(x,z)^rle Qiv) - Qiz) dz 

J d \g{z)\ 

and is stochastic if and only if Q(do) = oo. 

Let E — O = (0, oo) and a(x, 6) = x + 6. Then the operator P is the 
convolution operator with the measure v, i.e. if £i has density u and 9i has 
distribution v then Pu is the density of Ci + ^i- Assume that v has a density 
h. Then 

p(x,y) = l( QyX) (y)h(x -y), 
so that our evolution equation is 

du(t,x) dg(x)u(t,x) 



dt dx 



(p(x)u(t,x)+ h(x -y)ip(y)u(t,y)dy. 



Example 7. Friedman et al. [6] have considered stochastic aspects of gene 
expression following from bursts of protein production, and their formu- 
lations are special cases of our results. Identifying their w and c with 
w(x — y)dy = h(x — y)dy — 5 x (dy) andip(x) = k\c{x) and taking g{x) = —jx 
then our equation becomes identical, in a steady state situation to their 
Equation 6 

dx = / ^ y "> u *^ w ^ x ~ y "> dy - 
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Then, following [6], let h(y) = \e.~ y ^ b be the exponential distribution and 
let, in their first model, <p(x) = k±. Then the equation, as they have shown, 
has as a solution the density of the gamma distribution. In considering the 
second model of [6] that treated transcription factor regulation of its own 
transcription, we are let to consider the following function 

^ = kl TTx^ + kie ' 

As they have shown the corresponding density is given by 

a/a 



uJx) = cx a{1+e) - l e- x/b 



1 + x c 



where c is a normalizing factor. 
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